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r^P ' Abstract. Group classification of the three-dimensional equations describing flows of fluids 

with internal inertia, where the potential function W = W(p,p), is presented. The given 
equations include such models as the non-linear one-velocity model of a bubbly fluid with 
incompressible liquid phase at small volume concentration of gas bubbles, and the dispersive 
shallow water model. These models are obtained for special types of the function W(p, p). 
i-C . Group classification separates out the function W(p, p) at 15 different cases. Another part 

of the manuscript is devoted to one class of partially invariant solutions. This solution is 
f~| . constructed on the base of all rotations. In the gas dynamics such class of solutions is called 

the Ovsyannikov vortex. Group classification of the system of equations for invariant func- 
tions is obtained. Complete analysis of invariant solutions for the special type of a potential 
function is given. 

Key words: equivalence Lie group; admitted Lie group; optimal system of subalgebras; 
invariant and partially invariant solutions 
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1 1 Introduction 

OO ' 

The article focuses on group classification of a class of dispersive models [I 

^ \ p + pdiv(it) = 0, pu + Vp = 0, 

)W w _ (dW d (d 
bp ^ \ dp dt \ dp J \ dp 

where t is time, V is the gradient operator with respect to the space variables, p is the fluid 
density, u is the velocity field, W(p, p) is a given potential, "dot" denotes the material time 
derivative: f = % = ft + u ^f, and ^W- denotes the variational derivative of W with respect to p 
at a fixed value of u. These models include the non-linear one-velocity model of a bubbly fluid 
(with incompressible liquid phase) at small volume concentration of gas bubbles (Iordanski [2] , 
Kogarko [3], Wijngaarden [3]), and the dispersive shallow water model (Green & Naghdi [5], 
Salmon [6]). For the Green-Naghdi model, the potential function is [1] 

W(p,p)=p(3gp-e 2 p 2 )/6, 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http: / /www.emis.de/journals / SIGMA /symm etry2007,html 
1 See also references therein. 
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where g is the gravity, s is the ratio of the vertical length scale to the horizontal length scale. 
For the Iordanski-Kogarko-Wijngaarden model, the potential function is [T] 



£20 is the internal energy of the gas phase, c\ and C2 are the mass concentrations of the liquid 
and gas phases, n is the number of bubbles per unit mass, pio and P20 are the physical densities 
of components. The quantities c±, C2, n and pio are assumed constant. 

One of the methods for studying of differential equations is group analysis [7]. Many appli- 
cations of group analysis to partial differential equations are collected in Group analysis 
beside construction of exact solutions provides a regular procedure for mathematical modeling 
by classifying differential equations with respect to arbitrary elements. An application of group 
analysis involves several steps. The first step is the group classification with respect to arbitrary 
elements. This paper considers group classification of equations ([1]) in the three-dimensional 
case, where the function W pp satisfies the condition W pp ^ 0. Notice that for W pp = or 
W(p, p) = pif(p) + ip(p), the momentum equation becomes 

u + i)" Px = o. 

Hence in the case W pp = 0, equations ([I]) are similar to the gas dynamics equations. This case 
has been completely studied [9] (see also [TU]). 

The one-dimensional case of equations (pQ) was studied in [11]. As in the case of the gas 
dynamics equations there are differences in the group classifications of one-dimensional and 
three-dimensional equations. 

Another part of this paper is devoted to a special vortex solution. This solution was in- 
troduced by L.V. Ovsyannikov [12j for ideal compressible and incompressible fluids. This is 
a partially invariant solution, generated by the Lie group of all rotations. L.V. Ovsyannikov 
called it a "singular vortex". It is related with a special choice of non-invariant function. He 
also gave complete analysis of the overdetermined system corresponding to this type of partially 
invariant solutions: all invariant functions satisfy the well-defined system of partial differential 
equations with two independent variables. The main features of the fluid flow, governed by the 
obtained solution, were pointed out in [12]. It was shown that trajectories of particles are flat 
curves in three-dimensional space. The position and orientation of the plane, which contains 
the trajectory, depends on the particle's initial location. Later particular solutions of the sys- 
tem of partial differential equations for invariant functions were studied in [13 ^ 114 1 H6j. For 
some other models, this type of partially invariant solutions was considered in [17^ 118]. Exact 
solutions in fluid dynamics generated by a rotation group are of great interest by virtue of their 
high symmetry. The classical spherically symmetric solutions is one of the particular cases of 
such solutions. 

In this manuscript a singular vortex of the mathematical model of fluids with internal inertia 
is studied. Complete group classification of the system of equations for invariant functions is 
given. All invariant solutions for this system are presented. 

2 Equivalence Lie group 

Since the function W depends on the derivatives of the dependent variables, for the sake of 
simplicity of finding the equivalence Lie group, new dependent variables are introduced: 



W(p,p) 



p(c 2 p2o£2o(p2o) ~ 2irnp w R 3 R 2 ) 



where 




U5 =P, <j>i = W, <p 2 = W p , <j) 3 = W p , 
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where U4 = p and X4 = t. An infinitesimal operator X e of the equivalence Lie group is sought 
for in the form |19j : 

where all coefficients C Uj an d C^ k (* = 1>2, j = 1,2,3,4,5, = 1,2,3) are functions of the 
variables^ x,, ttj and 0^. Hereafter a sum over repeated indices is implied. 

The coefficients of the prolonged operator are obtained by using the prolongation formulae: 

C a " = DfC - u aJ DI^ (i = 1, 2, 3, 4), 
D- = d Xi + u^id Ua + {pxiWpj. + PxiWp^dwp, 

where a = (ai, a 2 , a 3 , 04) and /? = (/?i,/?2) are multiindices (oj > 0, > 0), 
(ai, a 2 , 03, "4), j = ("i + <5ij, «2 + <%, «3 + ^3j,a4 + 

«(a 1> a a ,as,a 4 ) = Qx? dx% 2 dxf 3 dt a * ' ^ lA) = dp~^dp~^' 

The conditions that IF does not depend on £, Xj, iij (z = 1, 2, 3) give that 

Cf = 0, Q=0, CT = °> C = (i = 1,2,3,4, j = 1,2,3, fc = 4,5). 

With these relations the prolongation formulae for the coefficients £ W P become: 
c w Pti = pe c w _ w^DfC - W^BfC* (i = 1,2), 

where 

Df = d p + Wf3,idw , b% = dp + w^dw, ■ 

For constructing the determining equations and solving them, the symbolic computer program 
Reduce [20] was applied. Calculations yield the following basis of generators of the equivalence 
Lie group 

Xf = d xi , X| = d X2 , X| = d X3 , XJ = td Xl + d Ul , Xf = td X2 + d U2 , 

Xq = td X3 + d U3 , X% = u 2 d U2 - uid U2 + x 2 d Xl - x x d X2 , 

Xf = u 3 d Ul - uid^ + x 3 d Xl - xid X3 , XI = u 3 d U2 - u 2 d us + x 3 d X2 - x 2 d xs , 

xf = d t , xf 1 = td t + x i d !H , xf 2 = dw, xf 3 = P d w , xf 4 = pd w , 

x f 5 = pd p + pd p + Wd w , Xf e = Xid x% + Uid Ui - 2pd p . 

Here, only the essential part of the operators Xf is written. For example, the operator Xf± 
found result of the calculations, is 

td t + Xid Xi - pdp. 

The part — pdp is obtained from X^ using the prolongation formulae. The symmetry opera- 
tors Xj (1 < j < 10) are symmetries of the Galilean groupH, which are independent of a potential 
function W(p,p). The symmetries corresponding to the operators Xf, Xf, Xf are the space 
translation symmetries, Xf , Xf , X| are the Galilean boosts, Xf, Xf and Xf are the rotations 

2 In the classical approach [7J Chapter 2, Section 6.4] for an equivalence Lie group it is assumed £J, fc = £^ = 0. 
Discussion of the generalization of the classical approach is given in [191 Chapter 5, Section 2.1]. 

3 This group is admitted by many systems of partial differential equations applied in Newtonian continuum 
mechanics. See, for example, [TJ [8] and references therein. 
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e a p, 


w' 


= e a W; 




w 


= W. 



and Xf Q is the time translation symmetry. The operator Xf± corresponds to a scaling symmet- 
ry, which is also admitted by the gas dynamics equations [7J. The symmetry corresponding 
to the operator Xf 6 applies for a gas with a special state equation [7]. Since the equivalence 
transformations corresponding to the operators Xf x , Xf 2 , ■ ■ ■ , Xf e are applied for simplifying 
the function W in the process of the classification, let us present these transformations. As the 
function W depends on p and p, only the transformations of these variables are presented: 

P 
P 
P 
P 
P 
P' 

Here a is the group parameter. 

3 Admitted Lie group of (CD) 

An admitted generator X of equations (pQ) is sought in the form 

x = t xi d xi + e 2 d X2 + e s 9 x , + eo t + c m d ui + c u2 q U2 + c 3 ^ 3 + c p a p , 

where the coefficients of the generator are functions of the variables x\, X2, £3, t, u±, 112, U3, p. 
Calculations showed that 

£ X1 = CQXit + C 4 t + C3X3 + X\C 7 + X1C1 + c 5 , 

£, X2 = C G X 2 t + Cut + X3C11 + X 2 C 7 + X 2 Ci - X\C\2 + C13, 

£Z3 = C6 x 3 t + C W t + C 7 X 3 + C1X3 - C n X 2 - C3X1 + C17, 

= c 6 t 2 + c 7 t + c 8 , ( p = (-3c 6 t + c 15 )p, 
( Ul = c 3 u 3 + c 2 u 2 - c 6 uit + c\U\ + c 6 xi + c 4 , 

C 2 Cll«3 - C G U 2 t + C\U2 - C 2 Ui + C 6 X 2 + C\2, 

C 3 = -ceu 3 t + citi 3 - cnu 2 - C3U1 + c 6 x 3 + ci 6 , 
where the constants q (i = 1, 2, . . . , 8, 11, 12, 13, 15) satisfy the conditions 
27c 6 p 3 (3W ppppP p + ^p PP p - 3W ppp p - W pp ) + 600W pp c 6 p 2 p 



+ 25p 3 (5W pppp p 2 (ci 5 - c 7 ) + 5W pppp ppci 5 + 18W ppp pc 15 
+ W pppP (28c 15 - 33c 7 - 10ci) + 18W pp {c 15 - 2c 7 - 2ci)) = 0, 
Wppppipj - ci 5 ) - ci 5/ 9Wppp + (2ci - C15 + 2c 7 )Wp P + 3c 6 W ppp p = 0, 
9Wppp P pp 3 ci5 + 40VFppp P/ 9 4 (c7 - ci 5 ) + Wppppp 3 /3(9c 7 - 49ci 5 ) - 9Wtfppf?f?CT 

+ 8Wpppp 3 (10 Cl - 17c 15 + 22c 7 ) + 2^pppp 2 p(9 Cl - 37c 15 + 9c 7 ) - 9W pppP 3 c 15 
+ 9W ppp pp 2 (c 15 - 2ci) + 56W>pp 2 (2ci - c 15 + 2c 7 ) + 9^ppp 2 (2ci - c 15 ) = 0, 
c 6 (5W ppp p + 3W>ppp + bW pp ) = 0. 



(2) 
(3) 



(4) 
(5) 



The determining equations ([2])-([5]) define the kernel of admitted Lie algebras and its exten- 
sions. The kernel of admitted Lie algebras consists of the generators 



YL = d xl , Y 2 = d X2 , y 3 = a 
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Y 4 = td Xl + d Ul , Y 5 = td X2 + d U2 , Y 6 = td X3 + d U3 , 
Y 7 = x 2 d X3 - x 3 d X2 + u 2 d U3 - u 3 d U2 , 
Y 8 = x 3 d Xl - xid X3 + u 3 d Ul - uid U3 , 
Y 9 = xid X2 - x 2 d Xl + uid U2 - u 2 d Ul . 

Extensions of the kernel depend on the value of the function W(p, p). They can only be operators 
of the form 

CiX 1 + CqXq + C7X7 + C15X14, 

where 

Xi = Xid Xi + Uid Ui , X 6 = t(td t + Xid Xi - Uid Ui - 3pd p ) + Xid Ui 
X 7 = Xid Xi + tdt, X 9 = x 2 d X2 + u 2 d U2 , X u = pd p . 

Relations between the constants c\, cq, C7, C15 depend on the function W(p,p). 

3.1 Case c 6 ^ 

Let cq / 0, then equation ([5]) gives 

5W pppP + 3W pppP + 5W PP = 0. 

The general solution of this equation is W pp = p~ 5 ^ 3 g(pp~ 5 ^ 3 ), where the function g is an 
arbitrary function of integration. Substitution of W pp into equation ([3]) shows that the function 
g = 2qo is constant. Hence, 

W = q p 2 p~ 5/3 + MP)P + MP), 

where the functions <p 2 {p) and fi(p) are arbitrary. Substituting this potential function in the 
other equations ©-([I]), one obtains 

3M" + ^ = 0. (07 + 201)^ = 0, C15 = -3(ci + c 7 ). 
If ip*2 = 0, then the extension of the kernel of admitted Lie algebras is given by the generators 

X§, X\ — 3Xi4, Xj — 3X14. 
If (p 2 ' = C 2 p~ 3 7^ 0, then the extension of the kernel is given by the generators 

Xq, X\ — 2X7 + 3X14. 

3.2 Case c 6 = 

Let cq = 0, then equation becomes 

-ci 5 a + (ci + c 7 )b + c 7 c = 0, (6) 

where 

a = pW ppp + pW ppp + W pp , b = 2W PP , c = pW ppp . 

Further analysis of the determining equations (J2])-([3|) is similar to the group classification of the 
gas dynamics equations [7]. 

Let us analyze the vector space Span(V), where the set V consists of vectors (a, b, c) with p 
and p are changed. If the function W(p,p) is such that dim(Span(V r )) = 3, then equation (|6|) is 
only satisfied for 

ci = 0, c 7 = 0, C15 = 0, 

which does not give extensions of the kernel of admitted Lie algebras. Hence, one needs to study 
dim(Span(y)) < 2. 
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3.2.1 Case dim(Span( V)) = 2 

Let dim(Span(V)) = 2. There exists a constant vector (a, (3, 7) 7^ 0, which is orthogonal to the 
set V: 

aa + (3b + jc = 0. (7) 
This means that the function W(p,p) satisfies the equation 

(a + i)pWppp + apW P pp = -(a + 2/3) W pp . (8) 

The characteristic system of this equation is 

dp dp dWpp 

{a + l)p ~^p~ -{a + 2(3)W P p 

The general solution of equation ([5]) depends on the values of the constants a, (3 and 7. 

Case a = 0. Because of equation (|7|) and the condition W pp 7^ 0, one has 7^0. The general 
solution of equation (jSJ) is 

W^(p,p) = <ppf t , (9) 

where k = —2/3/7, an d ¥> is an arbitrary function of integration. Substitution of ((9]) into ([6]) 
leads to 

ci5w' " <P(P) (2ci - (fc + l)ci5 + (A + 2)c 7 ) = 0. (10) 

If C15 7^ 0, the dimension dim(Span(y)) = 1, which contradicts to the assumption. Hence, 
C15 = and from (|10p one obtains c\ = — (k + 2)07/2. The extension of the kernel in this case 
is given by the generator 

- V X X + 2X7, 

where p = k + 2. 

If (A; + 2)(fc + 1) 7^ 0, then integrating one finds 

W(p, p) = ¥>(/>)p p + <pi(fi)p + <P2(fi), 

where <pi(p) and <p2(p) are arbitrary functions. Substituting this function W into ©-(jU) one 
has ip'2 = 0. 

If k = -2, then 

W(p, p) = ip(p) ln(p) + p(pi (p) + ip 2 (p) , 

and = 0, similar to the previous case. 
If k = — 1, then 

W(p,p) = (f(p)pln(p) + pipi(p) + tp 2 (p), 

and also ip' 2 ' = 0. 
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Case a/0, The general solution of equation ([9]) is 

Wpp(p,p)=<p(pp k )p x , (11) 

where k = — (1 + 7/a), A = —(1 + 2(3/ a) and <p is an arbitrary function. Substitution of this 
function into ([6]) leads to 

kop'z + kx<p = 0, 

where 

z = pp k , k = c 7 - cis(k + 1), fci = 2ci - ci 5 (A + 1) + 2c 7 . 
Since dim(Span(V)) = 2, one obtains that ko = and &i = or 

C7 = Cl 5 (fc + l), ci = c 15 (p - l)/2, 
where p = A — 2fc. Integrating (jlip . one finds 

W(p, p) = p*V(pp*) + pMp) + V2(p). (12) 
Substitution of CE2D into ©-gj) gives 

pip% + (2/c - A + 2)^2 = 0. 
Solving this equation, one has 

<f% = C 2 /" 2 , 

where C2 is an arbitrary constant. The extension of the kernel is given by the generator 
{p-l)X 1 + 2(k + l)X 7 + 2X 14 . 

3.2.2 Case dim(Span(V)) = 1 

Let dim(Span(V)) = 1. There exists a constant vector (a,(3,k) / such that 
(a, b, c) = (a, /?, k)B 

with some function B(p,p) / 0. Because W pp 7^ 0, one has that (3 ^ 0. Hence, the function 
W(p,p) satisfies the equations 

pW ppp + pWppp + (1 - 2a)W P p = 0, - 2jW pp = 0. 

The general solution of the latter equation is 

Wp P (p,p) = ip{p)p k 

with arbitrary function <p(p). Substituting this solution into the first equation, one obtains 
pip' (p) + (I - 2a + k)ip(p) =0, a = a/(3. 

Thus, 

Wpp = -q p k p\ (13) 
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where A = —(1 — 2d + k), qo is an arbitrary constant. Since dim(Span(V")) = 1, then go 7^ 0, 
A and k are such that A 2 + k 2 / 0. 

Substituting (fT3|) into ©, it becomes 

-c 15 (k + A + 1) + c 7 (k + 2) + 2ci = 0. 

Integration of (|13p depends on the quantity of fc. 

If (A; + 2)(fe + 1) 7^ 0, then integrating (|13p . one obtains 

WW) = -?opV + pMp) + ¥>2(p), p(p - 1) ^ o, 
where p = k + 2. Substituting this W into equations ((2])-(j3|), one obtains 

c 1 = (c 15 (p + X-l)-c 7 p))/2, 
with the function <p2(p) satisfying the condition 

C15W2" + ^2(- c i5(p + A - 2) + c 7 p) = 0. 
If = C2P fL 7^ 0, the extension of the kernel is given by the generator 

{l-p)X 1 +2{X u + <t>X 7 ), 
where (ft = (/i + \ + p — 2) /p. If ip' 2 ' = 0, the extension is given by the generators 

P X 1 -2X 7 , (p + A - l)Xx + 2X U - 
If k = — 2, then integrating (|13|) . one obtains 

W(p, P) = ~QoP X ln(p) + Wi (p) + <p 2 (p) , Qo + 0. 
Substituting this into equations ©-([I]), we obtain 

ci = ci 5 (A- l)/2, 
and the condition 

ci 5 (M" - Pa (A + 2)) + goA(A - l)(ci B - c 7 )p x - 2 = 0. 
If A(A — 1) = and 922 is arbitrary, then the extension is given only by the generator 

X 7 . 

If A(A — 1) = and ip 2 ' = C2p x+2 , then the extension of the kernel consists of the generators 

(A - 1)X X + 2X 14 , X 7 . 
If A(A - 1) + and ip' 2 ' = C 2 p x+2 - f A(A - l)pp x ~ 2 , then the extension is 

(A-1)X 1 +2(X 14 + (m+1)X 7 ), 

where 07 = (p, + l)cis. 

If fe = —1, then integrating (|13f) . one obtains 

w(p,p) = -<?op a pMp) + M(p) + 
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and substituting it into equations ©-([D), we obtain 

ci = (ci 5 A - c 7 )/2, 
and the condition 

C15P<P2 + ^2(- c 15A + Ci5 + Cy) = 0. 

One needs to study two cases. If ip? / 0, then the extension is possible only for tp2 = C2P" 11 / 0, 
where p = —A + 1 + C7/C15. The extension of the kernel is given by the generator 

(1 - p)X x + 2(/i + A - 1)X 7 + 2X U . 

If (p'2 = 0, then the extension of the kernel consists of the generators 

X\ — IXi, X\^ -\- XX \. 

3.2.3 Case dim(Span(V)) = 

Let dim(Span(I/)) = 0. The vector (a,b,c) is constant: 

(a, b, c) = (a, P, k) 
with some constant values a, (3 and k. This leads to 

Wpp = -2q , 

where go 7^ is constant. Integrating this equation, one obtains 

W(p,p) = -q p 2 + pMp) +Mp)- (14) 
Substituting (fT4|) into equation (|2])-(j4]), we obtain 

ci = (cis - 2c 7 )/2, 
and the condition 

If 7^ 0, then = C2P~ fl , where p, = 2c 7 /C15. The extension of the kernel consists of the 
generator 

(1 - p)X x + 2X 14 + fiX 7 . 

If (p'2 = 0, then the extension of the kernel is given by the generators 

X\ + 2X14, Xi — X 7 . 

The result of group classification of equations (pQ) is summarized in Table [TJ The linear part 
with respect to p of the function W(p, p) is omitted. Notice also that the change t — > — t has to 
conserve the potential function PF, this leads to <p>\{p) = 0. 

Remark 1. The Green-Naghdi model belongs to the class M7 in Table Q] with A = 1, p = 2 
and p = 0. Invariant solutions of the one-dimensional Green-Naghdi model completely studied 
in [2T]. 

Remark 2. The one- velocity dissipation- free Iordanski-Kogarko-Wijngaarden model has an 
extension of the kernel of admitted Lie algebras only for a special internal energy of the gas phase 
(class M3 (p = 2) in Table [T]), which corresponds to a Chaplygin gas £20 (P20) = 71/ P20 + Jo, 
where 71 and 70 are constants. 
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Table 1. Group classification of equations (JXJ) . 







Extensions 


Remarks 


Mi 


— qop~ 5 / 3 p 2 + LDoip) 


Xr, Xi — IXj -\- 3X14 


m" = C 2 p~ 3 ± 


M 2 




\)1 J_ _L ' -t 7 1 J_ rr 


(fn = 


M 3 


(f(p)pP + <p 2 


-pXi + 2X 7 


<p'> = 


M 4 


ip(p) lnp + ip 2 


X 7 


f'i = 


M 5 


pip(p) lnp + ip 2 


Xi - 2X 7 


¥> 2 ' = 


M 6 


P P <p(i>p k ) + V2 


(p-l)X 1 + 2{X 7 (k + l) + X u ) 


^ = c 2 ^- 2 


M 7 


-q p X p p + ¥2 


(1 - + 2(X 14 + (t>X 7 ) 


^' = c 2P ^ + 0, 

p(p - 1) / 0, 

4>={p + \ + p-2)/p 


M 8 




pX 1 - 2X 7 , 

(p + X- 1)X 1 + 2X14 


<& = 0, 
p(p -i)/o 


M 9 


-q p X lnp + ip 2 


x 7 


(P2(p) arbitrary, 
A(A - 1) = " 


Mio 




(A - l)Xi + 2X 14 , 
X 7 


^ = C 2P X+2 , 
A(A - 1) = 


Mu 




(A - l)Xi + 2(Xi4 + (m + 1)X 7 ) 


& = c 2P x+2 

do \ / \ i\ A— 2 

-fA(A-l) w A / , 
A(A - l) 


M12 


-q p x plnp + ip 2 


(1 - ^x)Xi + 2(/i + A - 1)X 7 + 2X M 


^ 2 = C 2 p-^ ± 


M13 




Xi — 2X 7 , X14 + AX 7 




M M 


-qop 2 + (f2 


(1 - p)X x + 2X14 + pX 7 


W = C^p^ / 


Mi 5 




X\ + 2X14, Xi — x 7 


<4 = 



4 Special vortex 

In this section a special vortex solution is considered. With the spherical coordinates |!2j : 

x = r sin 8 cos ip, y = r sin 8 simp, z = rcos8, 
U = u sin 8 cos (p + v sin sin (p + w cos 0, 
U2 = u cos cos 93 + v cos sin <p — w sin 0, 
E/3 = — u sin <p + v cos </?, 

the generators X 7 , Xg, X9 are 

X 7 = — sin ipdg — cos <p cot 88^ + cos <p(sin 8)~ 1 (U2djj 3 — Usdu 2 ), 

Xg = — cos ipdg — sin <p> cot #(9^ + siny?(sin 8)~ l (U2djj. A — U^djj 2 )i Xg = d v . 
Introducing cylindrical coordinates (H,u>) into the two-dimensional space of vectors (C/2,^3) 

U2 = H cos to, Us = Hsinu, 
the first two generators become 

X 7 = — sin ipdg — cos ip cot 88^ + cos c^(sin 8)~ 1 d LU , 
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Xg = — cos ipdg — sin tp cot 9d^ + sin <^(sin 9) d u . 

The singular vortex solution [12] is defined by the representation 

U = U{t,r), H = H(t,r), p = p(t,r), u = u(t,r,0,ip). 

The function uj(t,r,9,(p) is "superfluous": it depends on all independent variables. If H = 0, 
then the tangent component of the velocity vector is equal to zero. This corresponds to the 
spherically symmetric flows. For a singular vortex, it is assumed that H ^ 0. 

In a manner similar to [12] one finds that for system ([I]), the invariant functions U(t,r), 
Hit, r) and p(t, r) have to satisfy the system of partial differential equations with the two 
independent variables t and r: 

r 2 D p + p(r 2 U) r = pah, D U + p~ l p r = r~ 3 a 2 , 
D h = r~ 2 a(h 2 + 1), D a = 0, 

p = p(W p - pW pi) - W pp D p) + Wp P - W, (15) 

where a = rH, Dq = dt + Ud r , and the function h(t, r) is introduced for convenience during the 
compatibility analysis. 

The equivalence Lie group of equations (|15p corresponds to the generators 

XI = d t , XI = pd w , XI = 2td t - Udu - 3pd p - 5pd p - 3Wd w , 

xi = pd p + P dp + wdw , xi = xd x + udu + 2ad a + 2Wd w . 

Calculations yield that the kernel of admitted Lie algebras consists of the generator 
X = d t , 

extensions of the kernel can only be operators of the form 
hXi + k 2 X 2 + k 3 X 3 + k 4 X±, 

where 

Xi = td t - Udu - ad a + pdp, X 2 = t(td t + rd r - Udu - 3pd p - 5pd p ) + rdu - 3pd p , 
X 3 = 2td t + rd r - Udu - 3pd p - 5pd p , X 4 = pd p + pd p . 

The constants ki (i = 1, 2, 3, 4) depend on the function W(p, p). These extensions are presented 
in Table El 

4.1 Steady-state special vortex 

Let us consider the invariant solution corresponding to the kernel {^o}- This type of solution 
for the gas dynamics equations was studied in [14]. The representation of the solution is 

p = p(r), U=U(r), h = h(r), a = a(r). 

Equations (fT5|) become 

Up' + p(r 2 U)' = pah, UU' + p"V = r^a 2 , 
Uh! = r~ 2 a(h 2 + 1), Uol = 0, 

p = p(W p - Up'Wpp - W P pU(Up')') + WpUp' -W, p = Up'. (16) 
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Table 2. Group classification of equations ([15)) . 





W(p,p) 


Extensions 


Remarks 


Mi 


-qoP 2 p~ 5/3 + Pp 5/3 


X 2} X 3 


qoP 


M 2 


—qop 2 p~ 5 / 3 


A9, Xl , A 1 } 


QO 7^ 


M 3 


<p(j>)p" 


* 1 + (2-p)X 3 


p(p - 1) ^ 


M4 


-(QoP + VMP) + ( P2{p) 


-^1 + ^3 


c£>2 arbitrary 


M 5 


<p(p)pln(p) 


2Xx+X 3 




M 6 


p x <p{pp k ) + ip 2 (p) 


2Xi - (A - 2)X 3 , X 4 - kX 3 


^' = C 2P X ' 2 


M 7 


-q p X p p + <p 2 (p) 


2(pX 1 + 2(2p + p(X- p))X 3 
+{2-\){2X 1 + {2-p)X 3 ) 


y'l = C 2P » 

pip -i)/o 


M 8 


-qop X p p 


-2Xi + (2 - A)X 3 , 
(p - 2)X 3 - 2X 6 


p(p -i)/o 


M 9 


-q p x ln(p) +ip 2 (p) 


Xt + X 3 


9^2 arbitrary 
A (A - 1) = 


M10 




X 3 + X 6 , 2X l + (A - 1)X 3 


<ti = C2P X ~ 2 








A(A - 1) = 






^i + x - 2 -^x 3 + x§ 

+ CA(A-1) 1^-3 + A 6 j 


^=/- 2 ((/ ln(p)+/?) 

A(A — 1) f: 


M12 


A ■ 1 / ■ \ 1 / \ 


2Ai + Aa 3 

+(A-^-l)(X 3 + 2X 6 ) 


// /"~y // / n 

V? 2 = Cb/tr / 


M13 


-g p A pln(/)) 


2Xi + X 3 , X4 




M14 


-q p 2 +ip 2 {p) 


2Xi + x 3 - M x 6 


^' = C 2 p^ / 


M 15 


-qop 2 


Xi, X4 





In p3] it is shown that for the gas dynamics equations all dependent variables can be represented 
through the function h(r), which satisfies a first-order ordinary differential equation. Here also 
all dependent variables can be defined through the function h(r), but the equation for h(r) is 
a fourth-order ordinary differential equation. In fact, since H 7^ 0, from (|16p one obtains that 
U 7^ 0. Hence, a = otQ, where oto is constant. From the first and third equations of (|16p . one 
finds 

h! a (h 2 + l) 

p = Jtn ; =, (7 = — . 

In this case 



p = —aoRoh' 



y/W+T 



b! J 

and after substituting p and /) into the formula for the pressure, one has 
p = F(h,ti,h",ti"), 

where the function F is defined by the potential function W. Substituting representations 
of p, U and p into the second equation of (|15|) , one obtains the fourth-order ordinary differential 
equation for the function h{r). 
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4.2 Invariant solutions of ( TT5l) with W = — g p 2 P _5/3 + /3p 5/3 
System of equations (|15p with the potential function 

W = -q p 2 p" 5/3 + Pp 5/3 

admit the Lie group corresponding to the Lie algebra L 3 = {Xo,X2,X 3 }. 

If /3 = 0, then there is one more admitted generator X\. The four-dimensional Lie algebra 
with the generators {Xq,X\, X2,X 3 } is denoted by L4. 

The structural constants of the Lie algebra L4 are defined by the table of commutators: 





x 


Xx 


x 2 


X 3 


^0 





x 


x 3 











x 2 





x 2 









-2X 2 


x 3 












Solving the Lie equations for the automorphisms, one obtains: 



1 x 3 = x 3 + a x 2 , 1 ' 

, x 2 = x 2 + o 2 (xi + 2x 3 ) + a^xo, A 

' ^3 = ^3+ ^2^0, 



x 


= x e 


X-2 


= x 2 e 


Xo 


= x e 


X 2 


= x 2 e 



Construction of the optimal system of one-dimensional admitted subalgebras consists of using 
the automorphisms Ai (i = 0, 1, 2, 3) for simplifications of the coordinates (xo,xi,x 2 ,x 3 ) of the 
generator 

3 

X = J2 x i X r 

j=0 

Here k is the dimension of the Lie algebra (k = 3, 4). In the case L 3 one has to assume that 
the coordinate x\ = 0. 

Beside automorphisms for constructing optimal system of subalgebras one can use involutions. 
Equations fjl5[> posses the involutions E, corresponding to the change t — > —t. The involution E 
acts on the generator 

3 

X = J2 x i X r 

j=0 

by transforming the generator X into the generator X with the changed coordinates: 
E 

Here only the changed coordinates are presented. 



£0 = -x , 

X2 = -X2- 



4.3 One-dimensional subalgebras 

One can decompose the Lie algebra L4 as L4 = / © N, where / = L 3 is an ideal and N = {X±} 
is a subalgebra of L4. Classification of the subalgebra N = {X\} is simple: it consists of the 
subalgebras: 



Ni = {0}, N 2 = {Xi}. 
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According to the algorithm [22] for construction of an optimal system of one-dimensional 
subalgebras one has to consider two types of generators: (a) A = xoXq + X2X2 + X3A3, 
(b) X = X\ + xqXq + X2X2 + X3A3. Notice that case (a) corresponds to the Lie algebra L3. 
Hence, classifying the Lie algebra L4, one also obtains classification of the Lie algebra L3. 

4.3.1 Case (a) 

Assuming that xq 7^ 0, choosing 02 = — x^/xq, one maps X3 into zero. This means that X3 = 0. 
For simplicity of explanation, we write it as X3(^2) — ► 0. In this case ^2(^2) — ► X2 = X2 — x^ 2 /xq. 
If 5?2 7^ 0, then applying X2{A\) — > ±1, hence, the generator X becomes 

X2 + aXo, a = ±1. 

If X2 = 0, then one has the subalgebra: {Ao}. 

In the case xo = 0, if X3 / or X2 7^ 0, then, applying Aq, one can obtain xq 7^ 0, which leads 
to the previous case. Hence, without loss of generality one also assumes that X3 = 0, X2 = 0. 
Thus, the optimal system of one-dimensional subalgebras in case (a) consists of the subalgebras 

{A 2 ±A }, {X }. (17) 

This set of subalgebras also composes an optimal system of one-dimensional subalgebras of the 
algebra L3. 

4.3.2 Case (b) 

Assuming that xo 7^ 0, choosing 02 = —x^/xq, one maps X3 into zero. In this case X2(^2) - ► 
5?2 = X2 — £3(1 — xs)/xq. If x 2 7^ 0, then applying A\, and E2 (if necessary), one maps the 
generator X into 

X1+X2 + 7A0, 

where 7 7^ is an arbitrary constant. If x 2 = 0, then xo(^4o) ~~ > 0, and the generator X 
becomes X±. 

In the case Xq = 0, if 2x3 + 1 / or X2 7^ 0, then, applying Aq, one can obtain xo 7^ 0, which 
leads to the previous case. Hence, without loss of generality one also assumes that X3 = —1/2, 
X2 = 0, and the generator X becomes A3 — 2X\. 

Thus, the optimal system of one-dimensional subalgebras of the Lie algebra L4 consists of 
the subalgebras 

{A 2 ±A }, {A }, {X 1+ X 2 + 7X }, {X 3 -2X 1 }, {Ax}, 
where 7 7^ is an arbitrary constant. 

Remark 3. An optimal system of subalgebras for W = —qop^p^+ftp 5 with arbitrary f3 consists 
of the subalgebras (fTTl) . 

Remark 4. The subalgebra {X2 — Xq} is equivalent to the subalgebra: {A3}. 

4.4 Invariant solutions of X 1 + X 2 + jX 

The generator of this Lie group is 

A = 7X0 + X x + A 2 = (i 2 + t + i)dt + trd r - 3tpd p + (r - U(t + 1))% - ad a . 
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To find invariants, one needs to solve the equation 
XJ = 0, 

where J = J(t, r, p, U, a,h). A solution of this equation depends on the value of 7. 
Let 7 = p 2 + 1/4. In this case invariants of the Lie group are 

y = rs, V = s({(t + 1/2) 2 + p 2 )U -rt), R = ps~ 3 , A = ae^ arctan ^\ h, 

where 

s = ((t + 1/2) 2 + p 2 y 1/2 ^ arctan( ^ i) . 

The representation of an invariant solution is 

s (((t + 1/2) 2 + p 2 )U -rt) =V(y), p = s 3 R{y), a = Ae~^ ctM ^r\ h = h{y). 

Substituting the representation of a solution into (|15p . one obtains the system of four ordinary 
differential equations 

V = "I V + (Ah - 8Vy)/(Ay 2 ), h> = + A' = £, 

R'" = ( - ((8((3(4(44y + 5y)y - 19Ah)R + 308R'Vy 2 )R' 2 

- 3(88R'Vy 2 - 9AhR + 12(6V + y)Ry)R" R)Vq y - 9R 2/3 (4(4(2V + y)F 

- (4/i 2 + l)y V + (A - 4/iVy)A)li 3 )y - 18(8(R 2 / 3 Vy 3 + 4A/ig )V r y 

- (2/i 2 + l)A 2 g - 4(8(5F + y)V - (4// 2 + l)y 2 ) (? oy 2 ) J R , ^ 2 ))/(288 J R 2 y 2 go?/ 4 ). 
Let 7 = — /u 2 + 1/4. A representation of a solution is 

s(((t + 1/2) 2 - p 2 )U - rt) = V(y), a(t + 1/2 - p)^(t + 1/2 + p)~k = A(y), 
p(t + 1/2 - /U ) 3 ^ (t + 1/2 + p) 3a * = R(y), h = h(y), 

where 

y = rs, 8 = {t + l/2- p)'^(t + l/2 + p)'^, ai = ^f^, a 2 = ^-±1. 

4/x 4// 

In this case 

^ A/,-2^) /i , = A (^±l) A/= A 
R y 2 ' Fy 2 ' V 

R'" = {h28R"R'RV 2 q y i + 72R" R 2 Vq y 2 {-3Ah + 6Fy + y 2 ) - 6161? 3 V 2 q y 4 

+ 241? /2 l?ygoy 2 (19A/i - UVy - 5y 2 ) + 18R' R 2 {2R 2/3 V 2 y A - 8A 2 h 2 q 

- 4A 2 q + 32AhVq y - 40V 2 q y 2 - 8Vq y 3 - 4p 2 q y A + q y 4 ) 

+ 9R 2/3 R 3 y{4A 2 - AAhVy + 8V 2 y 2 + AVy 3 + 4p 2 y A - y 4 ))/{72R 2 V 2 q y 4 ). 

Let 7 = 1/4. A representation of an invariant solution is 

s((t + 1/2) 2 U - rt) = V(y), p = s 3 R(y), a = e 2 ^ 2t+ ^A(y), h = h(y), 

where 
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In this case 

a + (Ah-2Vy) tf _A(tf + I) A 
R y 2 V y 2 V' 

R'" = (528R"R'RV 2 q y 4 + 72R" R 2 Vq y 2 {-3Ah + 6Vy + y 2 ) - 616R' 3 V 2 q y 4 

+ 2AR /2 RVq y 2 (19Ah - UVy - 5y 2 ) + 18R' R 2 (2R 2/3 V 2 y 4 - 8A 2 h 2 q - AA 2 q 

+ 32AhVq y - 40V 2 q y 2 - 8Vq y 3 + q Q y A ) + 9R 2 / 3 R?y{AA 2 - AAhVy + 8V 2 y 2 

+ 4Vy 3 -y 4 ))/(72R 2 V 2 q y 4 ). 

These equations were obtained assuming that V 7^ 0. The case V = leads to 

A = 0, 2q R' - yR 5/3 = 0. 

4.5 Invariant solutions of X 3 — 2X\ 

Invariants of the generator 

X 3 - 2Xi = rd r - 3pd p + Udu + 2ad a 

are 

U = rV(y), p = r~ 3 R(y), a = r 2 A(y), h = h{y), 

where y = t. Substitution into equations (|15j) gives that the functions V(y), R(y), A(y) and h(y) 
have to satisfy the equations 

ti = A{h 2 + 1), A' = -2AV, R' = AhR, 

3{R 2 ' 3 + 6q )(V + V 2 ) = A 2 (4q (h 2 - 3) + 3{R 2 ' 3 + 6q )). 

4.6 Invariant solutions of Xi 

Invariants of the generator X\ 
X\ = tdt — Udu — 

are 

x, Ut, p, h, at. 
An invariant solution has the representation 

U = t- 1 V{y), p = R(y), a = t- 1 a{y), h = h(y), 
where y = x. Substituting into equations (fT5|) . one obtains 



R y 2 V y 2 ' V 

R'" = (132R"R'RV 2 q y 4 + 18R" R 2 Vq y 2 {-3ah + 6Vy + y 2 ) - l54R' 3 V 2 q y 4 

+ 6R ,2 RVq y 2 (l9ah - AAVy - 5y 2 ) + 9R' R 2 (R 2/3 V 2 y 4 - 4a 2 h 2 q - 2a 2 q 
+ 16ahVq y - 2W 2 q y 2 - Wq y 3 ) + 9R 2/3 R 3 y(a 2 - ahVy 
+ 2V V + Vy 3 ))/(l8R 2 V 2 q y 4 ). 

Here it is assumed that 7/0. The case V = only leads to the condition A = 0. 
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4.7 Invariant solutions of X 2 + X 

X 2 = t(td t + rd r - Udu - 3pd p ) + rdu. 
Invariants of the generator 

X 2 + X = (t 2 + l)d t + trd r - 3tpd p + (r - tU)d v 

are 

r^ + l)" 1 / 2 , U{t 2 + 1)V2 _ rt (t 2 + 1)-V2, + I) 3 / 2 , a, fc. 

An invariant solution has the representation 

£7^ + 1)V2 _ rt(t 2 + ^-1/2 = p = {t 2 + l)-3/2 E(y)) Q = ^ ^ = %) 

where y = r(t 2 + l) _1 /2. Substituting into equations (fl5|) . one has to study two cases: (a) V = 0, 
and (b) V ^0. 

Assuming V = 0, one obtains that A = 0, and the function R satisfies the equation 

2(5/?i? 4/3 - 9q )R' + 9yR 5/3 = 0. 
If V 7^ 0, then one obtains 

^ (M-2^) A(^+l) 
R y l V y 

R'" = {132R"R'RV 2 q y 4 + 54R" R 2 Vq y 2 (-Ah + 2Vy) - 154R /3 y 2 g 2/ 4 

+ 6i? /2 Wg y 2 (19A/i - 44Fy) - WR 1/3 R' R 3 (3y 4 + 9R'R 2 (R 2/3 V 2 y i - AA 2 h 2 q 
- 2A 2 q + l6AhVq y - 20V 2 q y 2 + 2q y A ) + 9R 2/3 R 3 y{A 2 - AhVy + 2V 2 y 2 
-y A ))/(18R 2 V 2 q y 4 ). 

4.8 Invariant solutions of X 2 — X 

Since the Lie algebra {X 2 — Xq} is equivalent to the Lie algebra with the generator {A3}, then 
for the sake of simplicity an invariant solution with respect to 

X 3 = 2td t + rd r - Udu - 3pd P 
is considered here. Invariants of the generator A3 are 

rt-Va, Ut 1 ' 2 , pi 3 ' 2 , h, a. 
An invariant solution has the representation 

U = t- 1 ' 2 V{y), p = t- 3 ' 2 R{y), a = a(y), h = h(y), 
where y = rt^ 1 / 2 . 

Substituting into equations (fT5j) . one has to study two cases: (a) V — y/2 = 0, and 
(b) V - y/2 ± 0. 

Assuming V — y/2 = 0, one obtains that A = 0, and the function R satisfies the equation 
2(20pR 4/3 + 9q )R' - 9yR 5/3 = 0. 
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If V — then one obtains 

V ( 10 in R ' . ^-(4V-3y)y A (h 2 + 1) 

V=(y/2-V)- + - 2 , h= {v _ y/2) y2 , A=0, 

= (2(y - 2V) 2 q y 3 (66R"R'Ry + 54R"i? 2 - 77i?' 3 y - 132R' 2 R) 
+ 9(y - 2V) 2 y 2 R 2 (R'R 2/3 y 2 - 20R'q + 2R 5/3 y) - 18R'R 2 y i q 
+ 6(y - 2V)ahy(18R"R 2 q y - 38R' 2 Rq y - 48R'R 2 q + 3R 2/3 R 3 y) 
- 72R'a 2 R 2 q {2h 2 + 1) - 40i? 10/3 Rl (3y 4 + 36R 2/3 a 2 R 3 y 
+ m u / 3 y 5 )/(18R 2 q y\y-2V) 2 ). 



5 Conclusion 

In this paper the complete group classification of the three-dimensional equations describing 
a motion of fluids with internal inertia (JT]) is given. The classification is considered with respect 
to the potential function W(p, p). Detailed study of one class of partially invariant solutions (the 
Ovsyannikov vortex) for a particular potential function is presented. This solution is essentially 
three-dimensional. 
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